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We present an exact diagonalization study of the self-energy of the two-dimensional Hubbard 
model. To increase the range of available cluster sizes we use a corrected t-J model to compute 
approximate Greens functions for the Hubbard model. This allows to obtain spectra for clusters 
with 18 and 20 sites. The self-energy has several 'bands' of poles with strong dispersion and extended 
incoherent continua with k-dependent intensity. We fit the self-energy by a minimal model and use 
this to extrapolate the cluster results to the infinite lattice. The resulting Fermi surface shows a 
transition from hole pockets in the underdoped regime to a large Fermi surface in the overdoped 
regime. We demonstrate that hole pockets can be completely consistent with the Luttinger theorem. 
Introduction of next-nearest neighbor hopping changes the self-energy stronlgy and the spectral 
function with nonvanishing next-nearest-neighbor hopping in the underdoped region is in good 
agreement with angle resolved photoelectron spectroscopy. 

PACS numbers: 71.10.Fd,74.72.-h,71.10.Ay 



I. INTRODUCTION 



Experiments on cuprate superconductors have shown a 
nontrivial evolution of their Fermi surface with hole dop- 
ing 5. In the overdoped compound Tl2Ba2Cu06-i-5 mag- 
netoresistance measurements[l[, angle- resolved photoe- 
mission spectroscopy (ARPES)[2] and quantum oscilla- 
tion experiments Q show a fairly conventional Fermi sur- 
face which is consistent with LDA band structure calcula- 
tions that take the Cu3d electrons as itinerant and which 
covers a fraction of the Brillouin zone of sa (1 — 5)/2. 
In the underdoped compounds the situation is more in- 
volved. ARPES shows 'Fermi arcs' |3| which however are 
probably just the intense part of hole pockets centered 
near (f,f). This is plausible because the sharp drop 
of the ARPES weight of the quasiparticle band upon 
crossing the noninteracting Fermi surface which must 
be invoked to reconcile the 'Fermi arcs' with the hole 
pocket scenario is actually well established in insulat- 
ing cuprates such as Sr2Cu2 02Cl2|5| and Ca2Cu02Cl2[l] 
where this phenomenon has been termed the 'remnant 
Fermi surface'. Moreover both the Drude weight in 
La2-2;Sr2;CuO40,ll| and YBa2Cu30y 8] as well as the in- 
verse low temperature Hall constant in La2-2;Sr2;Cu04'8:- 
[Tl| and YBa2Cu30j,[8] scale with 5 and the inferred band 
mass is constant troughout the underdoped regime and 
in fact even the antiferromagnetic phase 8j. This would 
be a exactly the behaviour expected for hole pockets. 
On the other hand, for 5 > 0.15 the Hall constant in 
La2-a;Sra;Cu04 changes rapidly, which suggests a change 
from hole pockets to a large Fermi surfaceQ. Quantum- 
oscillation experiments on underdoped YBa2Cu306.5fl2'- 
[Tsj and YBa2Cu40g[ll, show that the Fermi surface 
has a cross section that is comparable to <5/2 rather than 
(1 — 5)/2. Thereby the mere validity of the Fermi liq- 
uid description as evidenced by the quantum oscillations 
is clear evidence against the notion of 'Fermi arcs': the 
defining property of a Fermi liquid is the one-to-one cor- 



respondence of its low-lying states to those of a ficti- 
tious system of weakly interacting Fermionic quasipar- 
ticles and the Fermi surface of these quasiparticles is a 
constant energy contour of their dispersion and there- 
fore necessarily a closed curve in k-space. On the other 
hand the quantum oscillations cannot be viewed as ev- 
idence for hole pockets either in that both the Hall 
constant [ill and thermopower (Toj have a sign that would 
indicate electron pockets in the normal state induced 
by the high magnetic fields used in the quantum oscil- 
lation experiments. Thereby both, the Hall constant 
and the thermopower, show a strong temperature depen- 
dence and in fact a sign change as a function of temper- 
ature. At the same time neutron scattering experiments 
on detwinned YBa2Cu306.6 in the superconducting state 
show anisotropy in the spin excitations spectrum below 
30 meV and at low temperatures [20| . This indicates a 
rather complicated reconstruction to take place, possibly 
to a 'nematic' state with inequivalent x- and y-direction 
in the Cu02 plane. Such a nematicity which is also ap- 
parent in scanning tunneling microscopy [2lj must modify 
the Fermi surface in some way which may explain the un- 
expected sign. All in all the data may be interpreted as 
showing a change of the Fermi surface volume at around 
optimal doping from a small Fermi surface with a volume 
cx (5/2 to a large one with volume cx (1 — S)/2. 
Exact diagonalization studies of the t-J model have 
shown that the Fermi surface at hole dopings < 15% takes 
the form of hole pockets, 22 -24], that the quasiparticles 
have the character of strongly renormalized spin polarons 
throughout this doping range [25l - [27j and that the low en- 
ergy spectrum at these doping levels can be described as 
a Fermi liquid of spin 1 /2 quasiparticles corresponding to 
the doped holes [28j. A comparison of the dynamical spin 
and density correlation function at low[2l,|33| {S < 15%) 
and intermediate and high (5 = 30 — 50%) hole dop- 
ing moreover indicates [3 1| that around optimal doping a 
phase transition takes place. In the underdoped regime 
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spin and density correlation function differ strongly, with 
magnon-like spin excitations and extended incoherent 
continua in the density correlation function (HjH^I which 
can be explained quantitatively by a calculation in the 
spin-polaron formalism 32J . At higher doping, spin and 
density correlation function become more and more sim- 
ilar and both approach the self-convolution of the single- 
particle Green's function, whereby deviations from the 
self-convolution form can be explained as particle-hole 
excitations across a free electron- like Fermi surface fsT']. 
We thus expect a transition between a low-doping phase 
with a hole-pocket Fermi surface and quasiparticles which 
resemble the spin polarons realized at half-filling and a 
high doping phase with a free electron-like large Fermi 
surface. Here we want to further elucidate the issue of 
the Fermi surface and the possible transition between the 
large and small Fermi surfaces. To that end we study the 
electronic self-energy I](k, of the 2D Hubbard model 
by exact diagonalization. 



II. MODEL AND METHOD OF CALCULATION 

We study the Hubbard model on a two dimensional 
square lattice, defined by the Hamiltonian 
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4,aCj,a + H.C. 



^E^ 



(1) 



Here c\ ^ creates an electron with z-spin a in the orbital 
at lattice site i and denotes a summation over all 

nearest neighbor pairs. We set t — 1 and unless otherwise 
stated U/t = 10. The self-energy I](k, w) is defined by 
the Dyson equation 



where ek 
and 
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is the single particle Green's function at zero 
temperaturefs^. Here I^'q^'') and Eq^^ denote the 
ground state wave function and energy with N electrons. 
In the present study the Green's function for finite clus- 
ters is evaluated numerically by means of the Lanczos 
algorithm 1341 . 

Luttinger[33] has derived the following spectral represen- 
tation of the self-energy: 



S(k,w) 



.9k 



(4) 



In other words I](k, w) is the sum of a real constant gk 
(which is equal to the Hartree-Fock potential, see Ap- 
pendix A) and a sum of poles on the real axis. In the 



thermodynamical limit there may be both isolated poles 
and continua of poles Ck,<^- In ^ finite system, however, 
the poles Ck,;/ in principle always are discrete. Since the 
real part of E(k,ci;) assumes any value in [— cxi : oo] in 
between two successive poles, Ck.i^ and Ck,i/+i it follows 
that the equation 



(5) 



has exactly one solution in the interval [Ck.i/, Ck.i/+i]- If 
there is an energy interval with zero spectral weight - 
i.e. a gap - in the single particle spectral function it 
follows that there must be precisely one pole of the self- 
energy within in this gap. For example, the Hubbard-I 
approximation '35] for a nonmagnetic ground state corre- 
sponds to 



S(k,cj) = nU 



:j - (1 - n)U 



(6) 



where n is the density of electrons/spin. This is a single 
k-independent pole of strength cx at approximately 
the center of the Hubbard gap. In the neighborhood of a 
pole of I](k, w) the real part of the self-energy takes the 
form 



O'k.i/ 



Ck,i/ 



on the real axis, where Sr(a;) is slowly varying. If the 
residuum (Tk,!^ is large, the real part is large as well and 
no solution of a; — ek — 5R S(k, w) exists close to the 
pole. An isolated pole with large residuum thus 'pushes 
open' a gap of the spectral density in its neighborhood. 
On the other hand, if CTk.i/ is small the real part will de- 
viate from Tir{u}) only in the immediate neighborhood 
of Ck,i/- This implies that the corresponding solution 
of w — Ck = 5i S(k,aj) is pinned near Ck,i/. More- 
over, close to Ck,i/ the slope of the real part of the self- 
energy is large and negative, so that the spectral weight 
(1 — 9E(k, a;)/9a;)~^ is small. This rule will be seen fre- 
quently in the numerical spectra: an isolated pole with 
large residuum opens a gap in the single-particle spec- 
tral function around itself, a pole with small residuum 
has a pole of the single-particle Green's function with 
small weight in its immediate neighborhood. Finally, we 
note that a 'band' of poles of the self energy, i.e. Ck,i/ 
with V fixed, can never be crossed by a band of poles 
of the Green's function. Therefore, bands of isolated 
poles of S](k, w) define surfaces in the three-dimensional 
(fcj., fcj^, a;)-space which cannot be crossed by quasiparti- 
cle bands in the Green's function. The only exception 
would be a zero of the residuum, (Jk i^- 
As already mentioned we study the self-energy by com- 
puting the Green's function of finite clusters by means 
of the Lanczos algorithm. Thereby we encounter a tech- 
nical problem concerning the dimension of the Hilbert 
space. In a 4 X 4 cluster the dimension of the Hilbert 
space at half-filling (i.e. with 8 electrons of either spin 
direction in the cluster) is 65636900, in the half-filled 18- 
site cluster it is already 2363904400. Such large Hilbert 
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space dimensions make numerical calculations very diffi- 
cult. By contrast, the dimension of the Hilbert space of 
the Heisenberg model - which is equivalent to the Hub- 
bard model for large U /t - in the 4x4 cluster is only 
12870, in the 18-site cluster it is 48620. The Heisenberg 
model - and in the doped case the t-J model - thus are 
much easier to study numerically and in fact the largest 
cluster for which exact diagonalization studies for the 
Hubbard model have been performed [36l [37j is 4 x 4, 
whereas larger clusters are possible for the t-J model. For 
a study of the self energy, however, the t-J model cannot 
be used due to its 'projected' nature which for example 
implies the absence of the upper Hubbard band. 
On the other hand, various authors have derived effec- 
tive Hamiltonians which operate in the projected Hilbert 
space of the t-J model but reproduce physical quanti- 
ties of the Hubbard model to order t/C/jslSI. This 
is achieved by performing a canonical transformation 
which eliminates the part of the hopping term which cre- 
ates/annihilates double occupancies. The crucial point 
thereby is that not only the Hamiltonian itself, but all 
operators whose expectation values or correlation func- 
tions are to be calculated, have to be subject to this 
canonical transformation as well, which usually leads to 
correction terms of order t/U in all operators |4l| - [43| . If 
this is done consistently, however, very accurate approx- 
imate spectra for the Hubbard model can be calculated 
using 't-J-sized' Hilbert spaces which allows to treat more 
clusters and thus obtain additional information'43!]. So 
far this procedure has been performed only for the lower 
Hubbard band, because the study of the upper Hubbard 
band in the doped case requires a considerable number of 
additional terms in the Hamiltonian which describe the 
interaction between the doped holes and the double occu- 
pancy created in the inverse photoemission process (4]|. 
For the present study, however, the complete Hamilto- 
nian as given by Eskes et al.[AV\ has been implemented 
as computer code. A brief outline of the procedure and 
expressions for the corrected photoemission and inverse 
photoemission operators are given in Appendix B. This 
procedure allows to calculate approximate Green's func- 
tion for the Hubbard model over the entire doping range 
and on all clusters for which t-J model calculations are 
possible. For all these systems we evaluated the single 
particle Greens function ^ by the Lanczos method and 
obtained the self-energy from To illustrate the ac- 
curacy that can be expected Fig. [T] compares the single 
particle spectral function 

A(k,w) = - 5G(k,a;-i77) (7) 

TT 

for the true Hubbard model and the strong coupling 
Hamiltonian. While there are clearly some small differ- 
ences the strong coupling model reproduces the spectral 
function of the Hubbard model quite well. The devi- 
ations between the spectra calculated with the strong 
coupling model and the true Hubbard model are of order 
/U'^ for energies and /U'^ for the weights. This prop- 
erty in fact can be used to check the correctness of the 
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FIG. 1: Spectral density A{k,ijj) for a 4 x 4 cluster with 2 
holes computed with the true Hubbard model and with the 
strong coupling model. The ratio U /t = 10, r) = 0.1. 



strong-coupling-code be comparing energies and weights 
of peaks at diferent U jt. The agreement thus improves 
rapidly with decreasing tjU and already for U /t — 2Q the 
spectra become essentially indistinguishable. It is there- 
fore a useful check whether certain features of the spectra 
are robust with decreasing t/U. 

Lastly we mention that for the sake of analysis and ex- 
trapolation to the infinite lattice the energies ^k,!v and 
residua (Tk,iy of some poles will frequently be expanded in 
terms of tight-binding harmonics of k e.g. 
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70 (k) = 1 

7i(k) = 2 cos(fc3;) -I- 2 cos(fcj,) 

72 (k) ~ 4 cos{kx) cos{ky) 

73(k) = 2cos(2A:^) 2cos(2fcj/) (8) 



and similarly for the residua with coefficients aj. 
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FIG. 2: Spectral density A(k,ijj) and imaginary part of the 
self-energy S(k,a;) at half-filling. The figure combines spectra 
from the 16- and 18-site clusters, momenta from the 16-site 
cluster are marked by asterisks. The spectra are computed 
with 77 = 0.1, T,(k,uj) is multiplied by 1/50. Dots indicate the 
position oi = jj. + 2t{cos{kx) + cos(fcy)) for the respective 
momentum. 



III. RESULTS FOR THE GREEN'S FUNCTION 
AND SELF-ENERGY 



Figures [2] and [3] show the single particle spectral func- 
tion ^(k, a;) and the imaginary part i 3 E(k, w — if]) 
at half-filling. Particle-hole symmetry fixes the chemi- 
cal potential at /i = U/2. Figure [2] shows the entire en- 
ergy range of the lower and upper Hubbard band whereas 
Fig. [3] shows a closeup of the lower Hubbard band. Both 
Figures combine spectra from the 16 and 18 site clus- 
ter, which produces several k-points along each high- 
symmetry line. 

As expected S](k, w) shows an intense peak within the 
Hubbard gap which has a quite substantial dispersion. 
With the exception of the peaks at (tt, tt) and (0, 0) the 
dispersion of this peak is remarkably consistent with an 
inverted nearest neighbor dispersion, i.e. Ck — M = = 
2t{cos{kx) + cos(fcy)) which is indicated in Fig. [2l There 
are two possible interpreatations for the deviating be- 
haviour at (tTjTt) and (0,0). It may be that the domi- 
nant pole has a very rapid dispersion of its energy Ck.i/ in 



the neighborhood of these momenta. A second possibil- 
ity is that there are two poles with a smooth dispersion 
Ck,i/ but rapidly varying (Tk,!^- The peak at (7r,7r) and 
w « —6t would then belong to a second band of poles, 
which has appreciable residuum only near (tt, tt), whereas 
the residuum of the pole in the gap would suddenly drop 
to zero at {tt,tt). A calculation at U/t = 20 shows, how- 
ever, that the dispersion of the central pole including 
(tt, tt) and (0, 0) is almost exactly the same, which makes 
the interpretation in terms of a single central pole with 
rapid dispersion near these momenta more plausible. 
To understand the meaning of this form of the self-energy 
let us consider a self energy with a single dispersive pole 



S(k,..)^| 



V 



u 
2 



Ck 



(9) 



where the first term U/2 is the Hartree potential at half- 
filling. This yields the quasiparticle dispersion 



^k- 



(10) 



This is similar to spin density-wave mean-field theory 
which would be obtained by setting 77 = An ex- 

pansion of the type (jS]) gives the constant term ctq = 13.2 
for U/t = 10 and ctq = 84.1 for U/t = 20 so that we would 
obtain m = 0.53 for U/t = 10 and m = 0.84 for U/t = 20. 
On the other hand it is known from experiment that the 
dispersion of the spectral weight of the quasiparticle band 
at half-filling is not consistent with spin density- wave the- 
ory in that the spectral weight drops sharply in the outer 
part of the zonejl, Q . This is due to additional features 
in S(k, a;). 

The closeup of the lower Hubbard band in Fig. |3] re- 
veals additional structure in 'S(k.,uj). For most momenta 
there are two essentially dispersionless 'humps' at —At 
and —8.5t with a broad continuum in between them. 
The upper peak at —At shows some oscillation. It turns 
out, however, that the reason for this oscillation is that 
all peaks obtained in the 16-site cluster are shifted by 
~ —0.3t relative to those from the 18-site cluster. This 
shift and hence the entire oscillation may therefore be 
a finite-size effect. Together with the dominant peak 
within the Hubbard gap the peak - or group of peaks 
- at « — 4< encloses the quasiparticle band at the top of 
the photoemission spectrum. The intensity of the contin- 
uum is minimal at (0,0) and increases towards the zone 
boundary. Right at (7r,7r), however, the continuum is 
more or less absent. Also on this smaller scale I](k, w) 
thus shows a very rapid k-dependence in the neighbor- 
hood of (tt, tt). This behaviour is seen consistently in all 
clusters studied and is not an artefact of one specific clus- 
ter geometry. The 'band' of poles at ~ —At is the reason 
for the deviation from the simple spin-density wave form 
of the dispersion, Eq. (ITOl) . As can be seen in Fig. [3] the 
quasiparticle peak at (0, 0) is located immediately above 
the respective pole of the self-energy which has relatively 
small residuum. As discussed above, this implies a small 
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FIG. 3: Spectral density A(k,ijj) and imaginary part of the 
self-energy S(k,a;) at half-filling. The figure combines spectra 
from the 16- and 18-site clusters, momenta from the 16-site 
cluster are marked by asterisks. The spectra are computed 
with ri = 0.05, E(k, oj) is multiplied by 1/5 and at (7r,7r) by 
an extra factor of 0.25. 
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FIG. 4: Dispersion of the central peak of S(k,tj) in the Hub- 
bard gap for various doped systems. The lines are fits to 
tight binding harmonics with coefficients given in Table |T] 
The value U/t — 10 for all systems. 



the inverted nearest-neighbor hopping dispersion vanish 
with increasing U/t. 

The average residuum <to increases roughly as U^. Un- 
Hke the width of the dispersion, the weight of the central 
pole seems to be rather independent on doping. Figure H] 
moreover shows that already for two holes the peaks at 
(0,0) and (tTjTt) fit in smoothly into the dispersion and 
this holds true for all doped systems. 
In addition to this large peak for two holes a second band 
of less intense poles of E(k, w) appears in the neighbor- 
hood of (tTjTt). This can be seen in Fig. [S] which shows 
the spectral function and the self-energy for the lower 
Hubbard band. The pole in question starts out with the 
intense peak at w « /i — 2t at (tt, tt) and then rapidly 
disperses upwards. The respective peaks are very pro- 



weight of the quasiparticle peak itself. The dispersion- 
less band of poles at « —At thus reduces the bandwidth 

- because it cannot be crossed by the quasiparticle band 

- and also the spectral weight near (0, 0) and (tt, tt). 
We proceed to the hole doped case and consider the spec- 
tra for the cluster ground state with 2 holes, correspond- 
ing to S — 0.125 in 16 sites and 6 = 0.11 in 18 sites. The 
dominant pole within the gap still has a strong disper- 
sion although the bandwidth is reduced as compared to 
half-filling and the dispersion deviates from the simple 
inverted nearest-neighbor-hopping dispersion. This can 
be seen in Fig. U which shows the dispersion of this cen- 
tral peak for a few systems and in Table |T] which gives 
the corresponding parameters Q and ai for a variety of 
clusters. The following trends can be realized in Table HI 
for two holes the reduction of the bandwidth of the pole 
saturates at approximately 0.7 for large U/t. The devia- 
tions from the simple inverted nearest-neighbor-hopping 
dispersion seem to vanish in that limit. For 4 holes the 
same holds true, but the saturation value for the reduc- 
tion of the bandwidth is 0.4. The bandwidth of the cen- 
tral pole thus decreases with doping, the deviations from 
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-0.023 
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-0.420 
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80.216 


-0.312 


0.081 


0.094 
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-0.139 
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-0.438 



TABLE I: Expansion coefficients of the dispersion and 
residuum of the central pole in the gap for different systems. 
^: Data from the 20-site cluster. ": Data with t' /t = —0.4. 
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FIG. 5: Spectral density A(k,uj) and imaginary part of the 
self-energy E(k, oj) for the 16-site and 18-site cluster with 2 
holes. The spectra are computed with 77 — 0.05, E(k,a;) is 
multiplied by 1/5. 



FIG. 6: Spectral density A(k,uj) and imaginary part of the 
self-energy E(k, a;) for the 20-site cluster with 2 holes. The 
spectra are computed with -q = 0.05, E(k, is multiplied by 
1/5. 



nounced at (^,^) and (tt, ^), somewhat less clear at 
(tt, -1). This new band of poles can also be seen in Fig. 
[6] which shows the corresponding spectra for the 20-site 
cluster with 2 holes. At (tt, tt) itself the large peak is 
again at approximately u! = — 2t and at the two mo- 
menta (■^, ^) and ^) near (tt, vr) the intense peak 
is present as well. Since this upward dispersing band 
of poles can be identified only near (tt, tt) the residuum 
CTk.iy of this pole must have a strong k-dependence and 
decrease rapidly with increasing distance from (tt, tt). At 
(tt, tt) itself there is now also a broad incoherent contin- 
uum and the large peak at « /j. — 2t seems to have 
merged with this continuum. As was the case for half- 
filling the intensity of the continuum increases from the 
center to the edge of the Brillouin zone. Accordingly, 
the remnant of the free-electron band can still be seen at 
(0, 0) - this is the broad hump at — 4.5i in both, Fig. [5] 
and [6]- but is damped out for all other momenta. There 
is actually one difference between the 20-site cluster and 
the 16- and 18-site cluster: one might assign a third band 
of poles at the top of the incoherent continuum at w ~ — t 
in Fig. [5]- this band is completely absent in the spectra 
for 20-site cluster in Fig. [HI Figure [7] shows the spectral 
function and self-energy for the 20-site cluster with two 



holes and U/t = 20. This is qualitatively the same as in 
Fig. |6]but the bandwidth of the upward dispersing peak 
is reduced by approximately a factor of 2. The dispersion 
of this peak thus obviously scales with J = At^ /U and 
this is confirmed by other systems. 

The presence of an upward-dispersing band of poles of 
S(k, uj) near (tt, tt) and /i would be of crucial importance 
for the Fermi surface topology. Since the quasiparticle 
band cannot cross a band of isolated poles of E(k, w) this 
would force the quasiparticle band to bend downward [i^- 
l46j and thereby cut off the low energy inverse photoemis- 
sion weight at momenta around (tt, vr) in the energy range 
0.5t — >■ 2t in Fig. [5j This is a very plausible interpreta- 
tion, because - as shown in Ref. [43 - this low energy 
inverse photoemission weight is not part of a quasiparticle 
band, but spin-polaron shake-off. The downward bend- 
ing of the quasiparticle band in turn would lead to a band 
maximum and hence a hole-pocket-like Fermi surface - 
as found by exact diagonalization of the t-J model ^2j- 
[23 | and various version of Cluster Dynamical Mean Field 
Theorv [43 - l46j . As will be seen in a moment, this up- 
ward dispersing band of poles in the self-energy is indeed 
a special feature of the underdoped regime. 
Figure [5] shows the spectral function and self-energy 
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FIG. 7: Same as Fig. |6]but with U/t = 20. 



for 4 holes, corresponding to 5 — 0.25 in 16 sites and 
d = 0.22 in 18 sites and reveals a profound change in the 
self-energy. More precisely, the upward dispersing band 
of poles around (tt, tt) which was present in the under- 
doped case now has disappeared. The small peaks close 
to fj. which can be seen at (^, (7r,0) and (^,0) are 
probably a kind of finite size effect: at these momenta 
the quasiparticle peak is split between photoemission and 
inverse photoemission and since there is always a finite- 
size gap between phtoemission and inverse photoemis- 
sion spectrum in a finite system this results in a two- 
peak structure in the Green's function. This two-peak 
structure of the Green's function in turn necessitates a 
pole of the self-energy in between. In an infinite system, 
however, there is no splitting of a peak at /i by a finite 
amount, so that the respective peak in E(k, w) would be 
absent. Without these small peaks, however, S(k,ti;) has 
no significant peak above « fi — 0.8t. There are in fact 
some stronger peaks at the top of the continuum, particu- 
larly so at (tt, tt) and also at ^) and (§,§■) but if one 
wanted to assign a band this would rather have a shal- 
low maximum at (tt, tt) and then disperse downwards as 
one moves towards either (0, 0) or (tt, 0). The upward dis- 
persing band of poles in S(k, w) around (tt, tt) seen in the 
underdoped clusters thus is definitely absent, which im- 
plies a 'connected' nearest neighbor hopping band which 
starts out at w « — t at (0, 0) and reaches a; sa /i -f 2i at 
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FIG. 8: Spectral density A(k,uj) and imaginary part of the 
self-energy E(k, cj) for the 16-site and 18-site cluster with 4 
holes. The spectra are calculated with -q = 0.05 and S(k,tj) 
is multiplied by 1/5. 



(tt, tt). This band will produce a large Fermi surface but 
with a band mass that is enhanced by a factor of « 2.7. 
The further development with doping then is not really 
interesting any more: the central peak persists until hole 
dopings of 50% and becomes increasingly dispersionless, 
the lower Hubbard band stays similar to Fig. [5] 



IV. EXTRAPOLATION TO THE INFINITE 
SYSTEM 



Our next objective is to extrapolate the cluster results 
to the infinite system. It has to be noted beforehand 
that the results may not be expected to be quantitatively 
correct - this will become apparent by varios numerical 
checks - but rather give a qualitative picture. This is 
simply a consequence of the inavoidable limitations due 
to the small cluster size. 
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We represent the self-energy by the foUowing ansatz 

3 



E(k,w) ^ nU + J2 



3 

^Crc,A(k) log 



■ (11) 



The first term is the Hartree-Fock potential, the second 
term describes - for the underdoped case - the three dom- 
inant poles: the 'central pole' in the gap {v = 1), the 
upward dispersing pole near Q = (tt, tt) {v — 2) and 
the pole at the top of the continuum {i> = 3). The 
last term is the contribution from the incoherent con- 
tinua which we model by a constant spectral density be- 
tween k-independent limits e„iMi and Smax but with a 
k-dependent intensity crc(k). There are two such con- 
tinua in the lower Hubbard band, one below fi and the 
other above ^, and a third one for the upper Hubbard 
band. Since pole number 2, the upward dispersing pole 
near Q, can be seen only for a few momenta we ter- 
minate the expansion ([5]) after the second term for this 
pole, i.e. ^2,2 and ^3^2 are taken to be zero from the 
beginning. This pole moreover shows a strong variation 
of its residuum, crk,2, which rapidly decreases with the 
distance from (tt, tt) so that the pole cannot be identi- 
fied anymore for more distant momenta. Accordingly we 
approximate (Tk 2 as 



O'k.2 



/(k) = 4 cos^ ^ 



'/(k)/a 



cos 



2 / ky 



(12) 



Finally, the amplitude of the incoherent continua is writ- 
ten as 



(Tc(k) = CTc.O + CTc,! 1 ± 



71 (k) 



(13) 



where the — (+) sign refers to continua in the lower (up- 
per) Hubbard band. For the overdoped regime we use the 
same ansatz (jlip but without the upward dispersing pole. 
The pole at the top of the continuum (j^=2) has a rapid 
variation of its residuum as well so we use the expression 
([T^ . The coefficients which describe the dispersion of 
the central pole ly — 1 are given in Table U the remaining 
coeffcients are listed in Table [Hi Figure |9] compares the 
fitted self-energy in the underdoped case with the clus- 
ter spectrum. Fig. [TU] shows the same comparison for 
the overdoped vase. The agreement is not perfect but 
the fitted self-energy reproduces the essential features. 
The assignment of 'bands' in the self-energy clearly in- 
volves some degree of arbitrariness. It should also be 
noted that the dispersion of the pole 1/ = 2 has no sig- 
nificance in those regions of k-space where its residuum 
is small. Those parts which have large residuum, how- 
ever, appear to be fitted roughly correct. The fact that 
the band 1/ = 2 crosses the chemical potential leads to 




FIG. 9: S(k, oj) for U/t = 10 and two holes in the 16 and 
18-site cluster compared to the fit (|lip . The spectra are com- 
puted with rj = 0.05 and E(k,a;) is multiplied by 1/5. The 
dots indicate the dispersion of the poles = 2 and u — 3. 



additional inaccuracies: in a small cluster there is always 
an artificial finite gap between the photoemission and 
inverse photoemission spectrum, because the respective 
electron densities differ by a finite amount. This artificial 
gap can be up to 0.5i in the clusters studied and neces- 
sarily affects the dispersion of any band which crosses fi. 
This is certainly one reason for the inaccuracies of the fit 
for the band — 2. It therefore has to be kept in mind 
that the fitted self energies may not be expected to be 
quantitatively correct - rather the purpose of the fit is to 
illlustrate the consequences of the form of the self-energy 
in a more qualitative fashion. 

Next, we use the model self-energies to obtain approx- 
imate single-particle-spectra for the infinite system. For 
the underdoped (overdoped) regime we choose the elec- 
tron density per spin to be n = (1 — S)/2 with S — 0.12 
((5 — 0.24). We fix the chemical potential /i by demanding 
that the integrated spectral weight up to jj, to be equal 
to n: 



-E 



^(k, a;)(iw 



1 



(14) 



It turns out that the [i values obtained in this way de- 
viate only slightly (deviation < 0.2t) from the chemical 
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FIG. 10: Self energy for U/t = 10 and four holes in the 16 
and 18-site cluster compared to the fit. 



potentials of the cluster spectra. Figure [TT] then shows 
the single-particle spectral density for 5 = 0.12. The 
upper Hubbard band has been omitted because we rep- 
resented the self-energy in this energy range only by a 
continuum and did not attempt to fit any fine struc- 
ture. We note first that the spectral density in Fig. [TT] 
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^min ^max ^c,0 ^ c,l 






-6.0 -1.5 0.36 0.32 
1.0 6.0 0.30 0.0 





TABLE II: Coefficients of the model self-energy for 5 ^ 12% 
(top) and J K, 24% (bottom). The coefficients for the pole 
V = \ are given in Table U The constant terms Co,i/ and the 
edges of the continua are relative to ^i. 



FIG. 11: Single particle spectral density yl(k, lj) computed 
with the fitted self-energy pif) for the underdoped case 5 = 
0.12, the spectra are calculated with r\ — 0.05. 



is in very good agreement with the spectral density ob- 
tained by Quantum Monte-Carlo (QMC) simulations of 
the underdoped Hubbard model, see e.g. Fig. 9 of Ref. 
[isj . The two-band structure of the valence band, the fiat 
high-intensity part around (tt, 0) and the apparent near- 
est neighbor hopping band in the energy range [—t : t] 
are completely consistent with QMC. The intensity of 
the upper band in the photoemission spectrum is low at 
(0, 0) and increases as the Fermi energy is approached, 
whereas the lower band at — 3t has a high intensity at 
(0, 0) and rapidly looses weight as is moves away from 
this momentum. This is in agreement with the QMC 
spectra as well. 

Figure [T^] shows a Brillouin zone map of the spectral 
weight at /i. The Fermi surface obviously takes the form 
of a 'hole ring' along the surface of the antiferromag- 
netic Brillouin zone. Since the free dispersion ek is de- 
generate along the line (tt, O) — > (0, tt) there is no hole 
pocket but a 'hole ring'. Calculations for a single hole 
in the t-J model usually give - for moderate J/t - a very 
small dispersion along this line and a shallow maximum 
at (§■ f )(49l45^. If this band is filled with holes this 
results in hole pockets centered at (§,f). The reason 
for the maximum is hole hopping along a spiral path as 
first discussed by Trugman[49|. The present calculation 
either misses this fine detail or it is not relevant in the 
doped system so that no maximum exists and the pockets 
are deformed into a ring. In any way, the Fermi surface 
clearly is 'small' in that it covers only a tiny fraction of 
the Brillouin zone. This is ultimately the consequence of 
the upward dispersing pole of S(k, w) near {it,tt). The 
fraction of the Brillouin zone covered by the ring is 17.3%. 
This is much larger than the value of 5/2 = 6%. The 
latter value would be obtained if the doped holes are 




FIG. 12: Single particle spectral density at /i computed with 
the fitted self-energy (|ll|l for the underdoped case, S = 0.12. 
The value of rj is 0.05. 



FIG. 13: Single particle spectral density at /i computed with 
the fitted self-energy for the overdoped case. 



band. We set 



modelled by spin- 1/2 Fermions as suggested by exact 
diagonalization[28] and as predicted in a recent theory 
for lightly Mott insulators [53 • It is quite obvious, how- 
ever, that small changes in the parameters characterizing 
the fitted self energy may change this value strongly. The 
too large area of the ring thus probably simply shows the 
limited accuracy of the fit. One notable feature is the 
small spectral weight of the Fermi surface facing (tt, tt) - 
this is very similar to the 'remnant Fermi surface'. 
Next, we consider the overdoped case and set n — 0.38. 
Since there are no bands of poles close to and in par- 
ticular the band of poles near (tt, tt) is absent we expect 
a free-electron-like Fermi surface. This is indeed the case 
as can be seen from the Fermi surface map in Fig. 1131 
The fraction of the Brillouin zone covered by the large 
electron-like Fermi surface around (0, 0) is 0.382, which 
is in very good agreement with the Luttinger theorem 
when the carriers are electrons. Taken together the data 
thus indicate a phase transition in between underdoping 
and overdoping from a phase with hole pockets - or a 
'hole ring' in the present case - to one with a large Fermi 
surface. 



THE SELF-ENERGY IN THE 
UNDERDOPED REGIME 



We have seen that a key feature of the underdoped 
system is the presence of an additional upward dispersing 
band of poles in I](k, w) near (7r,7r), i.e. the band v — 2. 
In the following we discuss some consequences of this 



Ek + S(k,w + ^) = ek 



- Ck 



(15) 



where the first term on the r.h.s. is the sum of ek, ffk 
and the contribution from the other poles in E(k, w) and 
the second term represents the upward dispersing band 
I' = 2. We assume that ek is a smooth function of k and 
for simplicity neglect its frequency dependence. In the 
absence of the isolated pole ek therefore would be the 
quasiparticle dispersion. If additional poles of S(k, cj) 
are sufficiently far away we obtain the two bands 

wi,2 = ^ (?k + Ck ± V(ek-Ck)2+4(T) (16) 

which are shown in Fig. [TH There is a gap of 2y/a 
between these bands. Far from the crossing point Ck = Ck 
the bands take the form 



Wl,2 = Ek 



Ck - Ck 

Wl,2 = Ck - J- 

Ck - Ck 



(17) 



The two resulting bands thus partially trace the quasi- 
particle band ek and partially the dispersion of the pole, 
Ck- The spectral weight of the respective branches is 



1 



5ek ^ a 
dbJ (ck - Ck)^ 

^-1 = I I (^k-Ck)^ 

duj a 



(18) 



where the upper (lower) line refers to the band por- 
tion tracing Ck (Ck)- Since far from the crossing point 
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(ck — Ck)^ ^ cr the spectral weight Z assumes its usual 



value for the band portion tracing but is 



< 1 



for the band portion tracing This behaviour can be 
seen along (0,0) — )• (7r,7r) and along (7r,0) — > (tt, tt) in 
Fig. [11] The upward dispersing band of poles in I](k, w) 
also has a major significance for the validity of the Lut- 
tinger theorem i54i|. To see this we derive a slightly modi- 
fied version of the theorem which allows for an appealing 
physical interpretation. We consider 



S 



9S](k,cj) 



(19) 

where C is a curve in the complex w-plane which encloses 
the part of the real axis with a; < /i in counterclockwise 
fashion. All singularities of the integrand are located on 
the real axis. The first term in the second line will give 
the total electron number/spin. The second term has two 
kinds of singularties: poles and zeros of G(k, w). Near a 
pole we have 



1 



aG(k,a;) 



G(k,w) 9a; 
1 



(LJ-LJk,,y) 



(w - Wk,!/) 

whereas near a zero we have 

1 9G(k,a;) 
G(k,a;) 
1 



(20) 



Ck, 



CTkj(w - Ckj) 



+ 



- Ckj) 

It follows that 



(21) 



Ne = 2S' + 2^mk 

k 

mk = ^e(Ai-c^k,.) -5lQ(A^-Ckj). (22) 

If we assume that S" = we find that the number of 
electrons can be obtained by computing the number of 
'occupied' poles of the Green's function and subtracting 
the number of 'occupied' poles of the self-energy. A band 
of poles of the self-energy which crosses /i - such as the 
band v = 2 introduced in the fit of the self-energy in 
the preceeding section - therefore produces a 'negative 
volume Fermi surface' because the number of momenta 
within this surface has to be subtracted in the computa- 
tion of the electron number. 

To arrive at the known form of the theorem we note that 
- as discussed in section II - there is exactly one pole 
of G{k,uj) between any two successive poles of E(k, w). 
Moreover it is easy to see that there is always precisely 




(0,0) 



FIG. 14: Bare dispersion ek intersecting with a band of poles 
of the self-energy, C^k. 



one pole of G(k, lj) below the lowest pole of S(k, w). Ac- 
cordingly mk = if the topmost singularity below /i is 
a pole of the self-energy for the respective k-point and 
TOk = 1 if the topmost singularity is a pole of the Green's 
function. Since in the first case 3? G(k, /i) < whereas 
in the second case 3? G(k, /i) > we obtain 



= 2^e(5RG(k,A*)) 



(23) 



which is the 'generalized' Luttinger theorem given by 
Dzyaloshinskii[53|- The equivalence of (|22|) and (|23|) has 
previously been noted by OrtlofF et al. [56] . 

It is then easy to see that - contrary to widespread 
belief - hole pockets can in fact be completely consistent 
with the Luttinger theorem. We again consider Fig. 
1141 which shows a situation where the dispersion Ck is 
intersected by an upward dispersing band of poles of 
the self-energy, Ck, resulting in the two quasiparticle 
bands Wi^k and cj2,k- The lower of these bands, Wi^kj 
crosses the Fermi energy - indicated by the horizontal 
dashed line - and produces two Fermi level crossings 
at the points 1 and 2 which form the hole pocket. 
In between (0, 0) — > 1 we have rrik = 1 because the 
topmost pole below fi is the pole wi.k of the Greens 
function. Along 1 — > 2 mk = because no pole of either 
Green's function nor self-energy is below /i. Additional 
singularities at lower energies do not change this: since 
the topmost singularity below wi^k must be a pole of 
the self-energy, the total contribution to mk from all 
singularities including this one is zero. Along the short 
piece 2 — 3 we have again mk = 1 but along 3 — >■ (tt, tt) 
we have mk = because the topmost pole below fi 
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now is one of the self-energy, Ck- The piece 2 — >• 3 will 
be very short if the residuum a is small. The piece 
3 — )• (tt, tt) corresponds precisely to the 'negative volume 
Fermi surface' discussed above because here the topmost 
singularity is a pole of the self-energy. Assuming that 
the total volume of the hole pockets is Vbz (5/2 - as 
suggested by a recent theory of the lightly doped Mott 
insulator |53| - the fraction of the Brillouin zone outside 
the pockets is Vbz (1 — (5/2). Then, if the 'negative 
volume Fermi surface' is Vbz/2 the occupied part of 
the Brillouin zone in the sense of the Luttinger theorem 
would be Vbz(X ~ i-e. corresponding precisely 

to the electron density. Hole pockets therefore would 
be completely consistent with the Luttinger theorem 
if this is applied correctly. For example, evaluation of 
with the fitted self-energy for two holes - where 
the Fermi surface takes the form of a hole ring, see 
Fig. [12] - gives iVg = 1.032iV whereas the correct 
value would be iVe = 0.88iV. The deviation of « 15% 
shows the limited accuracy of the fitted self-energy, but 
is far smaller than the difference in Fermi surface volume. 



VI. COMPARISON TO ARPES EXPERIMENTS 

ARPES experiments on underdoped cuprate supercon- 
ductors show a number of interesting features and the 
next point is a comparison of the extrapolated cluster 
spectra to these experiments. Here an important point is 
to introduce longer range hopping terms. More precisely, 
we introduce an additional hopping term which connects 
next-nearest (i.e. (l,l)-like) neighbors. We choose the 
matrix element for this term to be t' = — 0.4t. This 
value is somewhat large but we simultaneously omit a 
hopping term between (2,0)-like neighbors because this 
would lead to a strong increase in the number of three-site 
terms in the strong coupling Hamiltonian. We performed 
the calculation only for the 16- and 18-site cluster with 
two holes because introduction of this additional hopping 
term increases the number of possible three-site combi- 
nations in the Hamiltonian considerably, so that the cal- 
culations become too difficult for the 20-site cluster with 
two holes and the 18-site cluster with 4 holes. It turns 
out that the ground states of the two clusters with the 
t'-term in the Hamiltonian have somewhat unusual quan- 
tum numbers: the ground state of two holes in the 16-site 
cluster has momentum (tt, 0) and spin S" = 0, the ground 
state of the 18-site cluster with two holes has momentum 
(^,0) and spin S* = 1. This means that the GS of the 
16-site cluster is twofold, that of the 18-site cluster 12- 
fold degenerate. This presents no real problem in that 
the expression (jSj should be viewed as the zero temper- 
ature limit of a grand canonical average, so that in the 
case of GS-degeneracy one simply has to average over all 
degenerate ground states. 

From Fig. |3]and Table U it can be seen that the disper- 
sion of the central pole within the gap (y = 1) is changed 
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FIG. 15: Spectral density A(k,tj) and imaginary part of the 
self-energy E(k, cj) for the 16-site and 18-site cluster with 2 
holes. The value t' = —OAt. The spectra are calculated with 
rj = 0.05 and E(k, cj) is multiplied by 1/5. 



somewhat by the presence of t' . In particular the devia- 
tions from the inverted nearest-neighbor-hopping disper- 
sion become stronger. Figure [TS] shows the single particle 
spectrum and self-energy in the presence of the i'-term. 
By comparison with Fig. [5] it is obvious that the t'-term 
introduces several pronounced changes in E(k, w): The 
residuum of the pole near (tt, tt) (z/ = 2) has decreased 
considerably, in fact this pole cannot seen anymore but 
takes the form of broad humps at the top of the inco- 
herent continua, whose upper edge has moved closer to 
fi. Instead there are now poles with large residuum at 
(tTjO) and (tt, ^). And finally there is no more pole at 
the top of the incoherent continuum and the intensity of 
the incoherent continua themselves has increased. To fit 
the self-energy we use the same ansatz PT|) but drop the 
pole v = 3 ed, the top of the incoherent continuum. The 
second difference concerns the strong pole near (tt, 0). 
We assume that this pole actually belongs to the band 
v = 2 and to model this we change the residuum of this 
band of poles by adding a second term: 

7d(k) = cos(fc^) - cos(fcj^). (24) 
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The dispersion of this pole is again expanded with re- 
spect to only two harmonics, the constant 70 (k) and the 
nearest-neighbor- hopping harmonic 71 (k). As already 
mentioned the poles near (tt, tt) cannot really be resolved 
in the calculated self-energies. Our main justification for 
keeping this band is the behaviour seen for t' = 0. 

Table Hill then gives the respective coefficients and Fig. 
1131 compares the numerical self-energy and the fit. The 
coefficient of 71 (k) is positive. If only the term cx ad were 
kept in ((M)) the exponential e^dl'^i replaced by unity and 
^o,iy be set to zero in this term, the resulting self-energy 
would be identical to the phenomenological self-energy 
introduced by Yang e< aL[57|. On the other hand the 
residuum clearly has an extremely strong k-dependence 
so that the exponential cannot be neglected in the present 
fit. 

Figure [TT] then shows the single particle spectral density 
obtained with the fitted self-energy. Along (0, 0) — >■ (tt, tt) 
the quasiparticle band disperses towards the Fermi en- 
ergy. Consistent with experiment, the intensity of the 
band thereby increases as [i, is approached. After cross- 
ing /i the band turns downward sharply and immedi- 
ately crosses \i again whereby its spectral weight drops. 
This is precisely the situation shown in Fig. [TH Along 
(0, 0) (tt, 0) the band disperses upward as well but 
does not reach /i. The ARPES spectrum thus shows a 
'pseudogap' but this is a trivial consequence of the Fermi 
surface being a hole pocket centered at (f , f ) (see be- 
low). Roughly at (^,0) there is a maximum of the dis- 
persion with high spectral weight and the band turns 
downward and looses weight beyond this point. This be- 
haviour may actually have been observed by Chuang et 
aZ. ^8] who interpreted this as indicating a Fermi level 
crossing at ~ (^,0). Chuang eA al. observed this be- 
haviour in an underdoped compound - see Fig. 2i of Ref. 
[ssj - but also in overdoped samples where it is unclear 
if it can be compared to the present calculation. In con- 
trast the spectral weight around (tt, 0) is small. This is 
consistent with experiment where a quasiparticle band 
around (tt, 0) is usually not observed in the normal state 
of underdoped cuprate superconductors. Along the line 
(tt, 0) — ?> (tt, tt) the band seems to disperse upward at 
first, but then again bends sharply and disperses away 
from 11. At the turning point the spectral weight drops. 
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FIG. 16: Imaginary part of the self-energy S(k, oj) for the 
16-site and 18-site cluster with 2 holes and Fit. The value 
t' = -OAt. 
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TABLE III: Coefficients of the model self-energy for U/t = 10, 
S ~ 10%, t' /t = —0.4. The constant term (0,2 and the edges 
of the continua are relative to fi. A comparison of the model 
self-energy and the numerical spectrum can be seen in Fig. 

m 



FIG. 17: Single particle spectral density A(k,uu) for 5 — 0.12, 
U/t — 10 and t' /t = —0.4 obtained by using the interpolated 
self-energy for the infinite lattice. 
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VII. DISCUSSION 




(0,0) (71,0) 



FIG. 18: Single particle spectral density at ^ for 5 = 0.12, 
U/t = 10 and t' /t = —0.4 obtained by using the interpolated 
self-energy for the infinite lattice. 



This is again due to the avoided crossing of the quasi- 
particle band with the upward dispersing band of poles 
of S(k, w) near (tTjTt), i.e. the band v = 2, m Fig. [HI 
A very similar behaviour has in fact been observed ex- 
perimentally in underdoped La2-2;Sr^Cu04, see Fig. 5 of 
Ref. |59|] . In this compound the quasiparticle band is suf- 
ficiently far from /x at (tt, 0) so that the absence of a Fermi 
level crossing is obvious for x = 0.05 and x = 0.10. A 
very avoided crossing has been observed as 'backbending' 
of bands in Bi2201 in Ref. [60]. If such an avoided cross- 
ing would occur sufficiently close to fi, however, it may 
look very similar to a true Fermi level crossing. Aparent 
experimental Fermi level crossings along (tt, 0) — t- (tt, tt) 
thus should be considered with care. 
The actual Fermi surface of the underdoped system is 
shown in Fig. [TH] and takes the form of a hole pocket, 
centered near (f , f ). The pocket is shifted slightly to- 
wards (0,0) and the part facing {tt,tt) has smaller spec- 
tral weight and less curvature than the part facing (0,0). 
The pockets covers 1.86% of the total Brillouin zone 
which would correspond - assuming twofold spin degen- 
eracy and four equaivalent pockets - to a quasiparticle 
density of 14.9%. This is close to the hole concentration 
of 12%. On the other hand the electron density as com- 
puted from the Luttinger theorem p3| is 1.25 so that the 
inaccuracies of the self energy clearly are substabtial and 
the close agreement for the quasiparticle density may be 
fortuitious. 



In summary we have presented an exact diagonaliza- 
tion study of the self-energy E(k, lj) in the 2D Hubbard 
model. Larger clusters than usual could be used be- 
cause instead of studying the true Hubbard model we 
considered its strong-coupling limit which requires much 
smaller Hilbert spaces. For dopings less than 30% - i.e. 
the doping region in which cuprate superconductivity 
takes place - several distinct features can be identified: 
first, a pole with large residuum c>c (U/t)^ and a disper- 
sion of width cx t in the center of the Hubbard gap, which 
is present throughout this doping range. Second, a pole 
with smaller residuum and an upward dispersion cx J 
around (tt, tt) which is present only in the underdoped 
regime. And third, several broad incoherent continua of 
width c>c t. The top of the lowest of these incoherent 
continua may be formed by a third pole. All features in 
the spectral representation of I](k, oj) show a pronounced 
k-dependence, both with respect to their dispersion and 
their residuum. This implies that in real space the self- 
energy is long-ranged and oscillatory. 
The key difference between the underdoped and over- 
doped system is the presence of a dispersive pole of 
I](k, w) around (tTjTt). This pole cuts through the quasi- 
particle dispersion and changes the Fermi surface topol- 
ogy completely. In the underdoped hole concentration 
range the Fermi surface takes the form of a 'hole ring' for 
t' — or hole pockets centered near (^, ^) for t' = —OAt 
and changes to a large free-electron-like Fermi surface in 
the overdoped case. We have shown that the hole pockets 
can be completely consistent with the Luttinger theorem. 
The residuum of the upward dispersing pole near (tt, tt) 
thus plays the role of an order parameter for the phase 
transition between the two different ground states. 
The single-particle spectra obtained with the fitted self- 
energies agree very well with both Quantum Monte Carlo 
simulations (for t' = 0) and ARPES on underdoped 
cuprates (for t' — OAt). In particular the spectra re- 
produce a hole-pocket-like Fermi surface of similar shape 
and location as in experiment and also the character- 
istic strong asymmetry of the spectral weight of the 
parts of the pocket facing (0,0) and (7r,7r). Quite gen- 
erally, the spectra show that all parts of the quasiparti- 
cle band which deviate from the noninteracting electron 
band structure have a very small weight. The reason is 
that these band portions closely follow the dispersion of 
a band of poles of S](k, w) and 'borrow' their spectral 
weight from the quasiparticle band. The data suggest 
in particular that some Fermi level crossings observed in 
ARPES along the line (tt, 0) — >■ (tt, tt) actually may not 
be true Fermi level crossings but sharp bends at the inter- 
section of the quasiparticle dispersion and a band of poles 
in S(k, Lu). An example where such pseudo-crossings can 
be clearly recognized is La2-a;Sra;Cu04[5^. In this com- 
pound the shift the quasiparticle band is sufficiently far 
from /i near (tt, 0) so that the absence of a Fermi level 
crossing is clear - if the quasiparticle band is closer to fi, 



15 



however, there pseudo-crossing may well be mistaken for 
a real Fermi level crossing. 

An interesting question is the physical meaning of the up- 
wards dispersing pole near (tt, tt). One immediate conse- 
quence of this pole is that the part of the inverse photoe- 
mission spectrum belonging to the lower Hubbard band 
consists of two disconneted components. The first com- 
ponent is the unoccupied part of the quasiparticle band, 
which forms the cap of the hole-pockets around (f , f ), 
the second component is a disconnected part around 
(7r,7r). This two-component nature of the inverse ph- 
toemission spectrum was discussed in Ref. [4?} where 
it was shown that the disconnected component around 
(tt, tt) actually consists of spin-polaron shakeofF, i.e. spin 
excitations which are released when a hole dressed by an- 
tiferromagnetic spin fluctuations is filled by an electron. 
The disconnected nature of the low energy inverse pho- 
toemission spectrum thus is a quite natural consequence 
of the hole pockets and the spin-polaron nature of the 
quasiparticles. 

The transition then may be understood as follows: in a 
Mott insulator the electrons are localized and retain only 
their spin degrees of freedom. Upon doping most elec- 
trons are still tightly surrounded by other electrons as in 
the insulator and thus remain localized - the ground state 
wave function therefore should optimize the energy gain 
due to delocalization of the holes and since the holes are 
spin-1/2 Fermions this can be achieved by forming hole 
pockets with a volume oc S around the ground state mo- 
mentum of a single hole i.e. ^). This picture leads to 
a simple theory [53| in which the number of noninteract- 
ing particles is equal to the number of doped holes. 
When the density of holes becomes sufficiently large so 
that the electrons are sufficiently mobile a phase transi- 
tion occurs to a state which optimizes the kinetic energy 
of the electrons themselves and this means the forma- 
tion of a Fermi surface with a volume (x 1 — d. A rough 
estimate for the hole concentration where the transition 
occurs would he Sc ~ z^^ with z the coordination number 
because then each electron has one unoccupied neighbor 
on the average. The vanishing of the upward dispers- 
ing pole in the self-energy then corresponds to this phase 
transition. Experimentally it seems that the phase tran- 
sition between pockets and free-electron-like Fermi sur- 
face occurs right at optimal doping. This suggests that 
superconductivity is also related to this transition. 
Acknowledgement: K. S. acknowledges the JSPS Re- 
search Fellowships for Young Scientists. R. E. most 
gratefully acknowledges the kind hospitality at the Cen- 
ter for Frontier Science, Chiba University. This work 
was supported in part by a Grant-in-Aid for Scientific 
Research (Grant No. 22540363) From the Ministry of 
Education, Culture, Sports, Science and Technology of 
Japan. A part of the computations was carried out at the 
Research Center for Computational Science, Okazaki Re- 
search Facilities and the Institute for Solid State Physics, 
University of Tokyo. 



VIII. APPENDIX A 

Here we proove that the real constant in ^ is equal 
to the Hartree-Fock potential. Considering the limit of 
large w, expanding the two expressions for the Green's 
function, ^ and ([3]), in powers of w"^ thereby using (HJ 
one obtains by comparing the terms cx a;~^: 



ek + 5k - (*o| 



(25) 



where [...]+ denotes the anticommutator. Using a stan- 
dard Hamiltonian of the form 



H 



k,cr 

2 ^ ^k+q,k,k'-q.k''^k+a.CT'^k'-a.<T''^k'.a'Ck.<T 
k,k' ,q cT,a"' 

this gives 

5k = XI K^lT'.k'.kK'^a') - K^ir',k,k'("k',<T). (26) 

k\cr' k' 

where (...) denotes the ground state expectation value. 

IX. APPENDIX B 

The canonical transformation which reduces the full 
Hubbard Hamiltonian to the strong-coupling Hamilto- 
nian takes the form 



O' = e^Oe-^ 



(27) 



where O denotes an operator in the original Hilbert space 
and O' is the transformed operator. The antihermitean 
generator S is 



EE 1(4. 



(28) 



The strong-coupling Hamitonian or corrected t-J model 
is obtained by transforming the Hubbard Hamiltonian 
according to (l27l) and discarding terms of 2'"* or higher 
order in l/U. The complete - and somewhat lengthy - 
Hamiltonian in given in equation (14) of the paper by 
Eskes et al. [4l| . The transformed version of the electron 
annihilation operator is 



e Cj^^e 



Eif( 4^-4.4) 



- (1 - 71.^4) (2ni,^ - l)dj-_^ ) 



(29) 
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where again terms of higher order in t/U have been ne- 
glected. By coUecting the terms which give a nonvanish- 
ing result when acting on a state without double occu- 
pancies and do not produce a double occupancy them- 
selves we obtain the operator for photoemission in the 
lower Hubbard band|4ll|: 

c»,t = c,,^ ~ 5Z If (^J.f^*'^ ~ S-.i-^r) • (30) 
j 

Since the transformed creation operator is the Hermitean 
conjugate of the operator for inverse photoemission 



with final states in the lower Hubbard band is just the 
Hermitean conjugate of pop. The operator for inverse 
photoemission with final states in the upper Hubbard 
band is obtained by taking the Hermitean conjugate of 
((29| and collecting terms which create a double occu- 
pancy but do not annihilate one: 

cit = ^«,t + X! "ff ( S"4 4i 4t 

j 

+ 4i + 4t (1 - "a)(2?^..t - 1) )■ (31) 
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